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Abstract The detection of the diffuse gas component of the cosmic web remains a formidable 
challenge. In this work we study synchrotron emission from the cosmic web with simulated 
SKAl observations, which can represent an fundamental probe of the warm-hot intergalactic 
medium. We investigate radio emission originated by relativistic electrons accelerated by shocks 
surrounding cosmic filaments, assuming diffusive shock acceleration and as a function of the 
(unknown) large-scale magnetic fields. The detection of the brightest parts of large (> lOMpc) 
filaments of the cosmic web should be within reach of the SKAl-LOW, if the magnetic field is 
at the level of a 10 percent equipartition with the thermal gas, corresponding to ~ 0.1 jj.G for 
the most massive filaments in simulations. In the course of a 2-years survey with SKAl-LOW, 
this will enable a first detection of the “tip of the iceberg” of the radio cosmic web, and allow 
for the use of the SKA as a powerful tool to study the origin of cosmic magnetism in large-scale 
structures. On the other hand, the SKAl-MID and SKAl-SUR seem less suited for this science 
case at low redshift (z < 0.4), owing to the missing short baselines and the consequent lack of 
signal from the large-scale brightness fluctuations associated with the filaments. In this case only 
very long exposures 1000 hr) may enable the detection of ^ 1 — 2 filament for field of view in 
the SKAl-SUR PAF Bandl. 


Advancing Astrophysics with the Square Kilometre Array 
June 8-13, 2014 
Giardini Naxos, Italy 


* Speaker. 


© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 


http://pos.sissa.it/ 









Filaments with the SKA 


Franco Vazza 


1. Introduction 

The Large-Scale Structure in the Universe comprises a complex network of filaments connect¬ 
ing virialized structures and separated by voids. Mildly nonlinear structures such as sheets and 
filaments should contain a wealth of information on the emergence of cosmic structures as individ¬ 
ually recognizable objects (e.g. galaxy groups and galaxy clusters). The observational detection 
of the warm-hot intergalactic medium associated with filaments is still challenging, and only a few 
works claimed a detection in X-rays (e.g. Richter et al. 2008), or more recently with the Sunyaev 
Zeldovich effect (Planck Collaboration et al. 2013). In the radio waveband, a few cases of radio 
emission at the possible crossroad of cosmic filaments have been reported (Bagchi et al. 2002; 
Kronberg et al. 2007; Pizzo et al. 2008; Giovannini et al. 2010; Farnsworth et al. 2013; Giovannini 
et al. 2015). However, such sources are likely to trace merger shocks in a filamentary environment, 
rather than stationary accretion shocks surrounding filaments. Numerical simulations suggest that 
filaments are surrounded by nearly stationary accretion shocks (M > 10), where the baryonic gas 
is shock-heated for the first time (Ryu et al. 2003; Vazza et al. 2011), and the detection of these 
shocks would confirm a critical piece of the warm-hot intergalactic medium (WHIM) paradigm 
(e.g. Dave et al. 2001). The expectations on the magnetic fields are more uncertain and lie in the 
range of ~ InG —0.1/iG, depending on numerical resolution of the cosmological simulation and 
as well as on the assumed seeding scenario (e.g. Dolag et al. 2008; Ryu et al. 2008; Donnert et al. 
2009; Vazza et al. 2014a). Provided that diffusive shock acceleration (see below) and large enough 
magnetic fields are present in filaments, the SKA could make the first detection of the cosmic web 
in the radio window (Keshet et al. 2004; Brown 2011; Araya-Melo et al. 2012). 

2. Simulated emission 

In Vazza et al. (2014b) we produced large cosmological simulations with grid code ENZO 
(Bryan et al. 2014), designed to study the dynamic properties of cosmic rays (CR) accelerated by 
cosmological shock waves. Here we mostly focus on a very large filament in our simulated vol¬ 
ume, about ~ 20Mpc in length and connecting two galaxy clusters (Fig.l). Given the size of this 
object, it gives us a first gross estimate of the chances of detecting filaments with the SKAl. Our 
forecasts of radio emission include two mechanisms: a) primary electron emission from relativistic 
electrons accelerated at shocks, assuming efficient diffusive shock acceleration (DSA, same model 
as in Hoeft & Briiggen 2007, Eq. 32); b) secondary electron emission from the relativistic elec¬ 
trons continuously injected by hadronic collisions in the filament volume (Dolag & EnBlin 2000), 
which are a by-product of our simulations (Vazza et al. 2014b). In general, the contribution from 
secondary electrons is everywhere negligible, i.e. only a few percent of the primary emission. 

The two main unknowns here are the electron acceleration efficiency at shocks {^e) and the mag¬ 
netic field within the filament (B); a discussion on these assumptions is given in Sec.^ The emitted 
radio power per frequency, V, scales as ly oc S ■ n ■ (where n and T are 

the post-shock density/temperature, S is the shock surface, and 5 is the slope of the CR-electrons 
energy spectrum). We can expect cosmic filaments to be a sensitive probe of cosmic magnetic 
fields, since for shocks in filaments (where 5^2 and B <C Bcmb) the previous formula reduces to 
ly oc -B^. In the following we assume E,e = 10^^ and uniform values of B = 0.01 — 0.3/xG within 
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Figure 1: First panel; project gas pressure map for a massive simulated filaments . Other panels: mock radio 
observations with SKAl-LOW, SKAl-MID and SKAl-SUR. The top panels shows the result for z = 0.05 
and B = Q.ljj.G. The lower panels are taken for the SKAl-LOW-C configuration (2 years survey at 300 
MHz and shows the result for different magnetic held strengths (B — 0.3, 0.1 and 0.03/rG) or redshift 
{z = 0.05 or z = 0.5). The additional blue circles represent the held of view of each array (not shown 
if larger than the image), while the red circles show the virial radius of galaxy clusters in the image. The 
contours (logio[Jy/arcsec^]) start from 3 times the nominal surface brightness (confusion) threshold for each 
array. 


the filament volume, and compute the detectable emission by SKAl-LOW (at two frequencies), 
SKAl-MID and SKAl-SUR. For each array we considered different possible strategies, including 
long 1000 hr exposures, 2-years surveys and the adoption point source subtraction at high resolution 
(~ 1”), followed by tapering on point-source subtracted data (at ~ 10”) to have a reduced resolu¬ 
tion (Table 1), since spatial-scale smoothing of the point-source subtracted data could significantly 
increase the signal-to-noise of the diffuse emission, though the power of identifying individual 
shock-structures would be lost (Brown 2011). For the SKAl-LOW observations we considered 
two possible central frequencies (110 and 300 MHz) and computed the confusion noise assuming a 
spectral index of —0.7, at full resolution. In the case of SKAl-MID and SKAl-SUR, we produced 
forecasts for a 1000 hr exposure or a 2 years survey. In all images we started from the full emission 
map obtained from the cosmological simulation, and then selected the specific frequency, field of 
view (FOV), resolution and sensitivity expected for each observation strategy, based on the SKAl 
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Table 1: Parameters for the investigated SKAl arrays and imaging strategies, based on the SKAl Baseline 
Design and on the SKA Phase 1 Performance Memo. For the tapering cases, we considered a point source 
subtraction at ~ 1”, followed by tapering on point-source subtracted data at ^ 10”. The sensitivities are 
expressed both as a function of the beam (1) and of the sky area (2); in the case of SKAl-LOW the sensitivity 
is set by the confusion noise, as in SKA Phase 1 Performance Memo by assuming a spectral index of —O.7.. 
For the SKAl-SUR and SKAl-MID we have considered a 30 percent fractional bandwidth for a 2 years 
all-sky survey, while for the SKAl-LOW we have considered an all-sky survey down to confusion noise, 
estimated from the rms confusion brightness at 1.4 GHz (Condon, private communication), with a confusion 
noise assumed as above. 


array 

configuration/strategy 

frequency 

[MHz] 

beam size 
[arcsec] 

field of view 
[degrees^] 

min. baseline 
[m] 

sensitivity 1 
[/rJy/beam] 

sensitivity 2 
[pjy/arcsec^j 

SKA-LOW-A 

full.res+conf., survey 

110 

6 

8 

45 

5.5 

0.134 

SKA-LOW-B 

full.res+conflOOO hr 

300 

3 

8 

45 

0.5 

0.0490 

SKA-LOW-C 

full.res+conf., 2 years 

300 

10 

8 

45 

9 

0.0794 

SKA-MID-A 

lOOOhr 

1400 

0.5 

0.49 

15 

0.09 

0.317 

SKA-MID-B 

2 years 

1400 

0.5 

0.49 

15 

5.8 

20.475 

SKA-MID-C 

1000 hr+tapering 

1400 

10 

0.49 

15 

0.14 

0.001 

SKA-MID-D 

2 years+tapering 

1400 

10 

0.49 

15 

8.8 

0.077 

SKA-SUR-A 

lOOOhr 

1200 

1.0 

18 

15 

0.4 

0.353 

SKA-SUR-B 

2years 

1200 

1.0 

18 

15 

3.8 

3.353 

SKA-SUR-C 

1000hr+tapering 

1200 

10 

18 

15 

0.5 

0.004 

SKA-SUR-D 

2years+tapering 

1200 

10 

18 

15 

4.9 

0.043 


Baseline Design of Mareh 2013 and on the SKA Phase 1 Performanee Memo. In addition, we 
FFT-transformed all images and removed spatial frequeneies smaller than the minimum baseline 
of eaeh eonfiguration, in the Fourier domain, so that the struetures with a too large angular seale 
eompared to the ehosen array are removed. In all eases, we implieitly assume a perfeet removal of 
the Milky Way foreground. 

3. Results 

Figure 1 shows a representative set of moek observations with various array eonfigurations of 
the SKA Phase 1. The top panels, show the performanees of the arrays in “deep exposure” mode 
(1000 hr), for a redshift z = 0.05 and assuming B = 0.1 flG ^ The SKAl-LOW-B array should 
deteet at least a signiheant portion of the filament at 3a, and traee some of the brightest emission 
knots eonneeting the two elusters in the field, as well as fhe big filamenlary aeerefion soufhwards 
of fhe mosf massive elusfer, reaehing ouf fo ~ 5 Mpe from fhe virial elusfer volume. Instead, wifh 
fhe SKAl-MID and SKAl-SUR arrays fhe defeefion of fhe filamenl will be impossible even wifh 
a 1000 hr poinfed observation, beeause of fhe smaller sensifivify and fhe laek of shorf baselines 
(af teas! fhe eenfral frequeneies of 1.4 and 1.2 GHz assumed here). The sifuafion gels even worse 
by eonsidering fhe shallower sensifivify reaehed in fhe 2-year survey. However, even wifh fhe 
SKAl-LOW fhe extended emission from fhe filamenl eonneeling fhe Iwo elusters would be barely 

*We remark that in this test we are only concerned with the possible emission from the WHIM from the filament, 
and that this emission model is likely overestimating the radio emission within cluster (e.g. radio relic emission), because 
of the assumed constant electron acceleration efficiency and magnetic fields. 
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detectable in its entire extent only for a higher magnetic field, e.g. B = 0.3/iG. A non-detection 
of such structures at low redshift with the SKA 1-LOW, would place robust upper-limits on the 
value of the intergalactic magnetic fields, of fhe order of < 0.03/xG, giving us imporfanf cluse 
on fhe efficiency of magnetic field amplifications in fhe WHIM. Af larger redshifls (z >0.1), fhe 
cosmological dimming makes fhe defecfion impossible, and also fhe removal of fhe confribufion 
from single radio galaxies becomes challenging. 

In order fo beffer quanlify fhe defecfabilify of fhe simulafed filamenf in fhe differenf arrays and using 
differenf possible imaging sfrafegies, we selected a smaller region inside fhe main filamenf (see fhe 
map in Fig.^, well oufside fhe surrounding clusfers oufskirfs, in order nof fo be confaminafed by 
clusfer accrefion shocks. We computed fhe fofal area of fhe source region above fhe 3a surface 
brighfness level, as well as fhe ratio befween fhe defecfable flux and fhe fofal one from fhe same 
region, as shown in Fig.§, for B = 0.3/xG or B = 0.1/iG. 

If B = 0.1/xG fhe SKAl-LOW is fhe only one giving any chance of defecfion. Using fhe planned 
array configurafion for fhe Phase 1 of SKAl-LOW our forecasl is fhaf ~ 30 — 60 percenf of fhe fofal 
flux from fhe filamenf would be defecfable, if mosf of fhis emission came from fhe brighfesf and 
densesf gas knofs wifhin fhe filamenf rafher fhan from fhe diffuse WHIM. Wifh fhe deep exposure 
of SKAl-LOW af 300 MHz (case B) fhe percenfage of defecfable flux approaches fo ~ 60 percenf, 
white fhe detectable surface goes to ~ 5 percent of the total. 

Based on the same computations, we can esimate that an increase in sensitivity of a factor x3 — 5 
will enable a clear detection of the full emission from such filament, up to several Mpc from the 
cluster centre, as well as the detection of the brightest emission spots from the filament at z = 0.5. 
This will be feasible in the Phase 2 of the SKAl-LOW, since the planned x20 increase of the 
maximum angular resolution (down to ~ 0.15“) in the 50 — 350 MHz range will lead to a great 
decrease of the confusion noise. 

On the other hand, in the case of SKA-MID and SKA-SUR even in Phase 2 the lack of short 
baseline would make it impossible to detect the bulk of synchrotron emission from the largest 
scales in filaments at the central frequencies we studied here. 

4. Discussion 

The forecast of radio emission from the diluted gas of cosmic filaments strongly depends on 
the combined value of ^e'B^, both of which uncertain. 

• The electron acceleration efficiency. The qualitative similarity to strong shocks in supernova 
remnants suggests that for the acceleration efficiency of electrons at strong shock 10~^ 
can be assumed (e.g. Edmon et al. 2011). However, the extrapolation of the DSA to the 
plasma condition of filaments is uncertain, since a too weak degree magnetisation would 
inhibit the formation of the magnetic field irregularities responsible for efficient particle 
scattering. Recent hybrid simulations of collisionless shocks suggest that this mechanism 
can still operate in presence of very low magnetisation, owing to the fast development of 
micro-instabilities at parallel shocks (Caprioli & Spitkovsky 2014). In general, any detection 
of radio emission associated to filament shocks will inform us that DSA (or some modified 
version of the theory) can operate down to this low density regime (n < lO^^cm^^), carrying 
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B=0.1nG 


surface fraction 
flux fraction 



B=0.3nG 



Figure 2: Fraction of the detectable area (red) and flux (blue) from the simulated filament, for B = 0.1/rG 
(left) and B = 0.3/rG (right). For a “detection” we considered a signal > 3(7, considering the sensitivites 
of Table 2. The inset shows the projected gas density of our filament and the selected white area used to 
compute the fraction of detectable radio flux. 

an important theoretical information. 


• The magnetic field. The expectations of MHD cosmological simulations lie in the range of 
~ 0.001 — 0.1/iG in filaments, depending on the numerical resolution and on the seeding 
scenarios. In the case of the giant filament analysed in this work, a value of B = 0.1/iG 
correspond to a magnetic energy density at the level of ~ 10 percent of the gas energy, 
which seems very hardly achievable by small-scale amplification of primordial seed fields 
alone (Briiggen el al. 2005; Ryu el al. 2008; Vazza ef al. 2014a). However, many addilional 
asfrophysical seeding scenarios of magnetic fields mighf add lo Ihe cosmological seed field 
(e.g. Donnerf el al. 2009). Any defection of coherenl emission on > 1 — 10 Mpc scales in 
filamenls will rule oul Ihe scenario in which Ihe magnefisalion of filamenls is driven by low- 
redshiff seeding of magnefic fields by galactic acfivily, owing lo Ihe difficulfy of covering 
such large scales wilh oulflows or diffusion. Al fhe same lime, Ihe defection of synchroiron 
from filamenls af high redshifl will give us imporfanl clues on Ihe efficiency of lurbulenf 
amplificalion of magnefic fields over cosmological epochs. 

Some insighl can also be gained from Ihe fact that no spectacular large-scale filament has been 
detected so far by LOFAR Cycle 0 observations. The present non-detection can be used to adjust 
our forecasts for SKAl. Following this conservative approach, we simulated the LOFAR LBA 
and HBA observations considering « 8 /r7/arcsec^ 1.7 /r7/arcsec^), as achieved in the current 
LOFAR-LBA (HBA) observations of nearby clusters, and a 23 x 16 arcsec^ (6 x tiarcsec^) beam 
In order not to be detected with present LOFAR Cycle 0 observations, B < 0.1 IJ.G, and we there¬ 
fore assume this conservative value for the maximum magnetic field in gianl filamenls, provided 
lhal DSA is al work in Ihese environmenls. As we have seen, even in Ihis case Ihe SKAl-LOW can 

^https://www.as tron.nl/radio-observatory/astronomers/array-configurations/. 
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Figure 3: Mock observations ofa ps 5 degree area within the filament, (as in the white selection of Fig.2, 
assuming a distance of z = 0.02), using the SKA simulator; 2000 hr pointing (left) and 60 hours pointing 
(right) at 140 MFIz with the SKAl-LOW (see text for more details). 



Figure 4: Left panel; projected gas temperature of a large cosmological volume at z = 0.5. Right panel; 
detectable emission by SKAl-SUR-A, using a 1000 hr exposure. The emission coming from inside the virial 
radius of all simulated clusters has been blanked out. The rectangular selection highlight a giant filament 
whose detection appears within reach. 


Still detect a significant portion of the filament emission for z < 0.1. 


5. Conclusions 

The detection of large cosmic filaments in radio at low redshift seems to be within reach by 
the SKAl-LOW, provided that the massive filaments are magnetised at the level of B ~ 0.1 jJ-G 
(~ 1 — 10 percent of equipartition with the thermal gas energy) and that DSA is at work in the 
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tenous WHIM environment. The SKA 1-LOW is always favoured because of its better sensitivity 
and sampling of short baselines compared to SKAl-MID and SKAl-SUR. In particular, the SKAl- 
LOW should be able to detect a significant portion of filament emission at low redshifts (z < 0.1), 
while for a better imaging of their full extent or for their detection up to z = 0.5 (or for lower 
magnetic fields) its sensitivity should be improved by a factor ~ 3 — 5. This can be achieved during 
the Phase 2 of SKA-LOW, thanks to the planned increased in the maximum angular resolution, 
which will enable to beat the confusion noise. On the other hand, our tests show that for most of 
cosmic filaments in the nearby (z < 0.1) Universe a detection with SKAl-MID or SKAl-SUR is 
made impossible by the lack of short baselines, for which even Phase 2 will not be beneficial. 

This work represents just a first step towards the realistic modelling of radio imaging of the 
cosmic web with the SKAl and the SKA2. Our schematic procedure of removing the missing 
baselines in the the Fourier plane is a very crude estimate of the actual UV coverage, and the effect 
of noise, confusion, finite bandwith etc. are not yet properly accounted for. In order to fill this gap, 
we produced a first set of simulated observations of our models with the most updated version of 
the simulated SKA performance In Figure 3 we present single pointing images at 140 MHz of a 
5 degree field of view centred on the faint filament region at z = 0.02 selected in FigJ^ for the 
R = 0.1 nG case. For the left image we used a 2000 hours of integration to achieve « 220 nJy/beam 
with 30 percent fractional bandwidth, while for the right image we simulated a 60 hour integration, 
achieving 1.3 /iJy/beam (in both cases using a 300 arcsec beam). The trends we recover with 
the use of the sophisticated tool provided by the SKA simulator overall agree with our previous 
estimates: while in a shallow survey mode only the brightest density peaks inside the filament can 
be marginally detected, the deep integration is able to pinpoint also the faint elongated structure 
of the filament, thereby suggesting its large-scale orientation and the presence of accretion shocks 
along several Mpc. 

Future deep exposures will also enhance the statistical chance of detecting even more massive 
filaments in large areas of the sky, not only with the SKAl-LOW if the targets are distant enough 
to have a sufficient baseline sampling even at higher frequencies. Figure ^ shows a test where 
we used a large simulated volume of (ISO^Mpc^) (using 1024^ cells) and simulated the detectable 
emission using a long 1000 hr exposure with SKAl-SUR (SKAl-SUR-A in Tab.l) and assuming 
a redshift of z = 0.5. In this case, we attached a magnetic field to the simulation by assuming 
an average 5 percent equipartition with the thermal energy within each cell. The field of view in 
the image approximately corresponds to 3 fields of view of the SKAl survey Bandl (rs 18 square 
degrees). Most of the detectable radio emission by the SKAl-SUR will originate from the cluster 
outskirts, yet in a few cases also the ’’tip of the iceberg" of radio emission coming from very big 
filaments should be detectable, like in the case of the hot filament within the rectangle in Fig.^. 
Our tests suggest that the redshift range 0.4 < z < 0.5 is the best to have some chance of detection 
with the SKAl-SUR PAF Bandl, while at z < 0.4 and z > 0.5 the emission will be filtered-out 
and suppressed by cosmological dimming, respectively. Based on a crude extrapolation of our 
data, we expect a detection of ~ 1 — 2 filaments of such giant filaments for every FOVs of SKAl- 
SUR, using ~ 1000 hr of exposure. Deep exposures in polarisation should maximise the chance 
of detecting these rare structures, since they are expected to be highly polarised owing to their 

^http://astronomers. skatelescope.org/wp-content/uploads/2014/06/SKAl_Science 
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Strong Mach numbers. Polarised observations are not dynamical-range limited, and warrant that 
the thermal noise will be reached independently of the presence of strong source within the field 
of view, thereby increasing the chances of imaging such low-surface brightness structures in the 
cosmic web. 
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